Loss of quantum coherence due to non-stationary glass fluctuations 
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Low-temperature dynamics of insulating glasses is dominated by a macroscopic concentration of 
tunneling two-level systems (TTLS). The distribution of the switching/relaxation rates of TTLS 
is exponentially broad, which results in non-equilibrium state of the glass at arbitrarily long time- 
scales. Due to the electric dipolar nature, the switching TTLS generate fluctuating electromagnetic 
fields. We study the effect of the non-thermal slow fluctuators on the dephasing of a solid state 
qubit. We find that at low enough temperatures, non-stationary contribution can dominate the 
stationary (thermal) one, and discuss how this effect can be minimized. 

PACS numbers: 03.65.Yz, 85.25.Cp 



The main hurdle on the way to implementation of a 
practical quantum computer lies in the design of a quan- 
tum system that is well isolated from the detrimental 
influences of environment, but at the same time acces- 
sible enough to allow coherent manipulation necessary 
to perform quantum computation. These are typically 
conflicting requirements since any external device that 
performs manipulation on a qubit is also a likely source 
of noise, which leads to decoherence. For example, metal- 
lic gates, which are commonly used to control solid-state 
qubits such as a superconducting Cooper pair box Q or a 
donor spin in a semiconductor [2j , as well as in ion traps 
0, are a source of electro-magnetic noise. This noise 
is generated by the low-energy particle-hole excitations 
in the metal, which in some cases is manifested as the 
Johnson-Nyquist Ohmic noise 0. 

Decoherence can also be induced by the low energy de- 
grees of freedom in the insulator surrounding the qubit. 
Recently, a number of studies were performed to an- 
alyze the role of the two-state fluctuators located in 
the insulators on the qubit dephasing and relaxation 
H H II H H 113 Such fluctuators naturally lead to 
charge noise. It is usually assumed that the fluctua- 
tors are in a thermal equilibrium, which is maintained 
via their interaction with the external environment and 
among themselves. However, very often the insulators 
used in the solid state qubits - S13N4, Si02 or AI2O3 
- are amorphous. It is well known that in amorphous 
materials, a large number of two-level systems exist with 
arbitrarily long switching times. In particular, they lead 
to the logarithmic time dependence of the specific heat 
[Tlj . Some of these slow two-level systems are frozen in 
high-energy states with energies significantly exceeding 
the nominal thermal energy and when they eventually 
relax to the ground state they stay there indefinitely. In 
this Letter we study the effects of such over-cooled two- 
level systems on the qubit dephasing. We find that for 



a typical glass they can give a significant contribution 
to the total dephasing rate, up to 1 GHz right after the 
system cool-down. 

The most common type of two level-systems (TLS) 
in structural glasses are the tunneling TLS, or TTLS. 
Their switching rate 7 is determined by either over- 
the-barrier thermally-activated tunneling, or under-the- 
barrier quantum tunneling. Thus the tunneling rate is an 
exponential of a uniformly distributed parameter charac- 
terizing the barrier strength. This generically results in 
the probability distribution V (7) = [7 log(7 max /7min)]~ 1 - 
An ensemble of TTLS with the energy splittings less than 
the thermal energy, Ei < kriT, naturally leads to the 1/f 
noise typical of glasses. The main contribution of 
this noise is concentrated at low frequencies ftto < kgT; 
therefore, it primarily affects the qubit dephasing (not 
relaxation), and moreover can be treated as a classical 
noise |y, |9j- The number of thermally active two-level 
systems scales linearly with the temperature (assuming 
flat distribution V(E) at small level splittings Ei), which 
for very low temperature would seem to imply that de- 
phasing should disappear. Here we argue, however, that 
at low enough temperature, slow non-thermal two-level 
fluctuators will dominate over the contribution from the 
thermal ones. That is because the number of non-thermal 
fluctuators is proportional not to the temperature but to 
the relevant band width T of TTLS, which can be up to 
20 Kelvin 0. 

The electromagnetic noise generated by the glassy 
metastable fluctuators can affect a variety of qubit ar- 
chitectures. However, for concreteness, here we consider 
the effect of the electric field fluctuations generated by 
the TTLS switching on a superconducting qubits. Elec- 
tric field fluctuations couple to the qubit charge. In the 
qubit charge basis, the Hamiltonian in the presence of 
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the charge noise X(t) is 



H 



E c 



Ej 



X(t) 



(1) 



where a z and o~ x are the Pauli matrices acting on the 
space of the qubit states {0, 1} and Ec and Ej are the 
qubit charging and Josephson energies. Diagonalizing 
the noise-free part of the Hamiltonian and keeping only 
the longitudinal coupling to the noise (assuming that its 
contribution is negligible at the qubit frequency fkv 
E^^E, 

relaxation) the Hamiltonian becomes 



Ej and thus cannot effectively lead to qubit 



H 



E 



Ec 
2E 



X(t)a z 



(2) 



Note that TTLS can be also present in the insulat- 
ing parts of Josephson contacts. Their switching then 
leads to fluctuations of the Josephson energy Ej 0,0], 
the relevant contribution to the Hamiltonian (J2J) being 
(Ej /2E)X(t)a z . This contribution is important only 
near degeneracy point, Ec — 0, where it can be taken 
into account in a similar way. 

For slow fluctuators the noise can be treated classi- 
cally and it is a sum of individual contributions X(t) = 
J2 v i z i(t), where u, is the strength of the coupling be- 
tween the qubit and the i fluctuator, and Zi(t) describes 
the time evolution of the fluctuator. To single out the ef- 
fect of non-thermal fluctuators, we neglect here the ther- 
mally excited fluctuators, such that E < ksT. Without 
loss of generality, we set for the excited fluctuator 2 = 1 
and for the ground state z = 0. 

The dephasing effect of a particular realization of the 
fluctuators' dynamics on the qubit can be defined as 



n< 



if*dt'v iZi (t') 



(3) 



It describes the phase drift between the and 1 sates 
of the qubit due to the random switching events in the 
environment, as would be measured in a free induction 
decay (FID) experiment 0. Here v = (Ec/E)v. The 
dependence on the biasing point Ec / E is routinely used 
to disentangle the charge noise acting on qubit from the 
phase/flux noise. To determine the statistical dephasing, 
we need to average over the time evolutions z, (t) , the pa- 
rameters of the fluctuators (the switching rates 7), and 
the coupling strengths v. Assuming uncorrelated fluctu- 
ators, we can perform the evolution averaging indepen- 
dently for all Zj. For an exponential relaxation from the 
excited state with switching rate ji, V(l\l) — e~ 7i '*', the 
average single-fluctuator contribution to dephasing is 



dt'^e'^'e^' +e~ 7 ' t e ii ' t 



IVi - Ji 



(4) 



This expression corresponds to the average dephasing 
factor due to the i th fluctuator, which one would ob- 
tain in an experiment repeated with the identical initial 
conditions (fluctuator i in the excited state at t = 0). 
Accordingly, the dephasing factor due to an ensemble of 
uncorrelated fluctuators is 



*(*)=n^(*)- 



(5) 



When many fluctuators interact with the qubit simul- 
taneously, one can evaluate the dephasing factor by per- 
forming the average over the distribution of {Si,7t}- We 
make the following assumptions: (i) Fluctuators are dis- 
tributed in 3D according to the Poisson distribution, 
dN(r) = n4irr 2 dr, where n is the concentration; (ii) 
The coupling to qubit is of dipolar form v(r) — A/r 3 , 
where A can be either positive or negative; (iii) There 
is no correlation between 7 and v; (iv) the distribution 
function for the switching rate has the tunneling form 
p {l) = (£7)~\ where C = ln(7, nax /7 min ). Under these 
assumptions, the full distribution function (normalized 
to the total number of fluctuators N = nV) , 



V(v, 7 ) =N^-^e 



rj/\v\ 



n 



4irA 



(0) 



The physical meaning of 77 is the coupling between the 
qubit and the fluctuator located at the boundary of the 
volume V (the "weakest" fluctuator). 

To evaluate dephasing factor in Eq. J5J, we use the 
Holtsmark procedure |17j for averaging the logarithm of 
the product, 



/C* lu (i) = -ln^ lu (i)« / dvd 1 V(v, 1 )g^ u {v n \t), 

(7) 

where g FlB (v,j\t) = 1 - fo£+ u D (t) + V'-™(*)]/2 includes 
averaging over the sign of coupling v. According to our 
assumptions, the distribution function does not depend 
on the sign of v. Thus we can perform the sign average 
in ip and limit the integration over v to positive values 
only. From Eqs. J3J and ©, we obtain for t > 



J 7min ' 



-7t\ 



For j max t <C 1 this yields, 

K FID (t«7 m ax)=/^typi, 

while in the other regime 

K FID (7ma X « t « J m l) - V typ t hi ( — 

\7minl 



(8) 



(9) 



(10) 



The characteristic coupling strength that deter- 
mines the decay of coherence is defined here as 
v typ = (E c /E)2tt 2 AN '/3V£. Note that both short and 
long-time decay laws @ and (|10f> are nearly exponential. 
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The free induction decay is the simplest measure of de- 
phasing. However, even a fluctuator that never switches 
on the timescale of the experiment will give a contribu- 
tion to FID dephasing if we use definition J3J). Moreover, 
experimentally it is impractical to reset ( "re- initialize" ) 
the states of the fluctuators. This dependence on the fluc- 
tuator initial-state can be easily eliminated by measuring 
the dephasing relative to the initial qubit level splitting, 
which includes contributions from the two-level fluctua- 
tors. The result is the following definition for dephasing 



m 



which yields 



7n 



IC IC \t) = -v typ / -4 (l-jt-e-y*) , 



and in the short and long-time limits 

relA/A _ / ^typTmax^ j t <C 7max 



/C rel (i) 



ln(7maxt) Vty P t , 7m i x < t < 7 mi 1 n 



(11) 



(12) 



(13) 



With this definition, the long-time dephasing is similar 
to FID, however, in the short time limit the dephasing is 
quadratic in time and thus dramatically suppressed. 

Example with the relative dephasing shows that sub- 
traction of the systematic qubit frequency shift dramati- 
cally reduces dephasing. However, direct implementation 
of such subtraction protocol is not possible, except for the 
slowest fluctuators that remain frozen from one experi- 
mental run to another. Alternatively, a similar subtrac- 
tion effect can be obtained by means of spin echo 0] 
protocol. In the echo experiment, the static phase accu- 
mulation is eliminated by introducing n pulse on qubit 
at time t/2 and performing measurement at time t, 



i^S-S t,2 dt'x(t')-S t t/2 dt'x(t') 



Integration over v distribution now yields 



/C cch °(i) = i typ I % (l + e-f* - 2e~^ 2 ) (15) 



(14) 



and 



IC echo (t) 



max'' j ^ < ^«» 7max 

d typ t In 2, 7 m * x < t < 7"^ 



(16) 



Indeed, the result is qualitatively similar to the one for 
the protocol (|11|) . Notice also that the long-time dephas- 
ing for all protocols, including FID, is the same. 

In thermal equilibrium, the maximum switching rate 
7max(T) is determined by the interaction between a 
TTLS with electrons or phonons. The estimates have 
been obtained for amorphous metals jl9j | and dielectric 
glasses [2(J- At a given temperature T, these estimates 



can be are 7max(T) xs yTjh where \ = 0.01 — 0.3 
for amorphous metals, and 7max(T) ~ T 3 /hT 2 where 
T c = 15 — 30 K (depending on the elastic parameters) 
for dielectric glasses. The minimal switching rate, 7 m i n , 
is actually not known. Logarithmic heat release from 
structural glasses was observed during many hours, see 
e. g. Ref. |2l| and references therein. The low-frequency 
noise in disordered materials has 1/f spectrum down to 
any observable frequencies. 

The non-stationary case considered here, however, re- 
quires careful interpretation of 7 max and 7 m i n - The upper 
limit of the switching rate, 7 max , is determined by the 
time delay, t<i, since the preparation of the system (e. g., 
since the cool-down, or some other strong manipulation 
that can reset the fluctuators). Indeed, the fluctuators 
with 7 > 1/rd are likely to decay before the the start 
of the measurement. This naturally leads to a cut-off 
7max = l/^rf. The simplest assumption about the lower 
cut-off on relaxation rate is that it is the same as in the 
equilibrium case. However, it can also can also depend 
on history. For instance, if the system was warmed-up 
to temperature To for a period of time t p , then the fluc- 
tuators with energies less than To, and relaxation rates 
faster than l/r p will get "recharged." To include this pos- 
sibility, we define 7 m ; n = 1/t p . 

Since the measurement time in a typical experiment is 
less than Ta we can concentrate only on the "short-time" 
limit, t <C 7max- Then, 



K{t < r) 



Vty-pt , 

tityp 4-2 



FID 

echo . 



(17) 



It was shown recently that the short-time limit of the 
echo decay has the desirable self-averaging property, that 
is the ensemble average (calculated here) indeed corre- 
sponds to the typical result from the a single sample, 
averaged over repeated experimental runs [y, 12^ . On 
the other hand, FID and the long-time limit of echo do 
not self-average due to strong mesoscopic fluctuations in 
the positions of the nearest fluctuator (s)!^- 

We now crudely estimate the magnitude of dephas- 
ing assuming that the qubit environment is an insulat- 
ing glass, such as SiC>2 or A1 2 C>3. The concentration 
of TTLS's in such systems is about Vo = N jVL ~ 
10 32 — 10 33 cm _3 erg _1 [l3j. The energy distribution 
of TTLS is uniform, with the band width of about 
20 Kelvin. The band width is much larger than the 
typical temperature 10 — 50 mK in the qubit dephas- 
ing/relaxation experiments, and therefore only a small 
fraction of TTLS are thermally excited and contribute 
to the telegraph noise capable of dephasing the qubit. 
On the other hand, the fraction of metastable fluctua- 
tors that do not decay during the time since the cool- 
down is M = VoVT \n(T p /Td). Here T is the relevant 



bandwidth of the nonequilibrium fluctuators, e. 
temperature before final cool-down. 



the 
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The coefficient in the dipolar interaction for a TTLS 
that interacts with its mirror image on the surface of 
the qubit is A ~ e 2 d 2 , where the typical dipole mo- 
ment, ed, corresponds to one electron charge displaced 
by d ~ 1 A. Assuming the bandwidth of the frozen 
TTLS to be To ~ 1 K we obtain as a crude estimate, 
Utyp/fi ~ 2n 2 e 2 d 2 ToT a ~ 10 8 - 10 9 s -1 . Introducing now 
the dephasing time as the time during which the log of the 
average phase factor becomes equal to one, K-iT^) = 1, 
we find Tf ID ~ {E/E c )v£ ~ (E/E c ) • (10~ 8 - lO" 9 ) s 



and f 2 echo 



r d f 2 FID . 



It is instructive to compare these results with the 
contribution to dephasing from the thermally excited 
fluctuators. Let us make this comparison for the echo 
signal. The effective number of thermal fluctuators is 
less than number of the frozen ones by the factor ~ 
T/Tq <C 1. Since vt yp is proportional to concentration, 



frozen /..thermal 



"typ 



typ 



T /T > 1. For thermal fluctuators 



T: 



ccho =mm 1 U tyP 



thermal /..thermal 



typ 



■ 7max(T) 



At experimentally relevant temperatures, in dielectric 
glasses one can expect 7 max (7 1 ) <C vlyp Tmed . For example, 
at T = 50 mK, v t ^ rmal ~ 10 6 s" 1 and 7 max (T) - 10 3 s _1 . 
Therefore, 

^frozen / rp \ 1/2 

2 i i / T^n — 



^thermal ^ 



Assuming that 7 max (T) = T 3 /T 2 h we get 



frozen 



^thermal j; ^ ftfl. 



1/2 



(18) 



This relation shows that at 



r d « KT T 2 /T A 

frozen fluctuators provide more decoherence than ther- 
mal ones. Assuming T = 1 K, T c = 20 K, T = 50 mK, 
we conclude that the crossover occurs at ~ 10~ 2 s. 
While this delay time is much shorter that a typical time 
from a cool-down to an experiment, it is much longer 
than the typical experimental inverse repetition rate.jlj 
In case the manipulation procedure includes application 
of significant bias or gate voltages, e.g. at the measure- 
ment stage, this can lead to "charging-up" of the non- 
thermal fluctuators. In particular, a temporary swap- 
ping of ground and excited states of a fluctuator can be 
induced by the gate electric fields. 

In conclusion, we analyzed the effect of slow non- 
thermal metastable glass fluctuators on qubit dephas- 
ing. We found that at low enough temperatures, the 



non-thermal dephasing can exceed the contribution of 
the thermal fluctuators. This effect can be reduced by 
appropriately designing manipulation sequences in order 
to avoid re-exciting of the metastable fluctuators and by 
choosing sufficiently long delay time between consequent 
manipulations. 
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